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We study the nature of the disorder-induced quantized conductance, i.e., the phenomena of topological Ander- 
son insulator (TAI) induced in HgTe/CdTe semiconductor quantum well. The disorder effect in several different 
systems where anomalous Hall effect exist, is numerically studied using the tight-binding Hamiltonian. It is 
found that the TAI phenomena also occur in the modified Dirac model where the quadratic corrections k^a- 
is included and electron-hole symmetry is kept. It also occurs in the graphene system with the next nearest- 
neighbor coupling and staggered sublattice potential. Comparison between the localization lengths of the 2D 
ribbon and 2D cylinder clearly reveals the topological nature of this phenomena. Furthermore, analysis on the 
local current density in anomalous quantum Hall systems where the TAI phenomena can or can not arise reveals 
the nature of TAI phenomena: the bulk state is killed drastically and only the robust edge state survives in a 
moderate disorder When the edge state is robust enough to resist the strong disorder that can completely kills 
the bulk state, TAI phenomena arise. 

PACS numbers: 73.23.-b, 73.43.-f, 73.20.At, 72.15.Rn, 



I. INTRODUCTION 

It is known that the two-dimensional (2D) noninteracting 
system with the quadratic dispersion relation is an Anderson 
insulatori^ 

XXX 

In the presence of strong spin orbit coupling (SOC)24, 
the externali^ or internal^ magnetic field, the metallic state 
is present. Due to the crossing of the mobility edge^, a 
metal-insulator transition (MIT) occurs at the critical disorder 
strength Wc where the localization is divergent. Recently, Li et 
al^ find that the disorder drive either a metallic state or an ordi- 
nary insulating state to the topological insulator (it is so called 
topological Anderson insulator) in the HgTe/CdTe quantum 
well, which has been numerically confirmed by Jiang et el'-. 
Using the effective medium theory, the mechanism of the TAI 
is explained as the crossing of a band edge rather than a mo- 
bility edge by Groth et al.'" Although there are many inves- 
tigations focused on the TAI, there are still some unanswered 
questions. For instance, what leads to the band edge cross- 
ing, or, what's the nature of the TAI phenomena (i.e., disor- 
der induce quantized conductance)? Furthermore, except for 
HgTe/CdTe quantum well, is there any other systems that have 
TAI phenomena? What's the necessary condition to generate 
TAI phenomena? 

As we know, the edge states in topological insulator-ii are 
"helical" states, i.e., the direction of propagation is tied to the 
electron spin rather than electron charge, the edge states lo- 
cated in the opposite edges are protected by the time reversal 
symmetry and each contribute e^/Zi to the conductance. In 
principle, we can focus only on one edge state which is tied 
to spin up or spin down to study topological insulator. The 
sub-system related to each individual spin contribute to an in- 
dividual anomalous quantum Hall effect. In the following, we 
will deal with the sub-system related to the spin up in topolog- 
ical insulator and treat it as an individual anomalous quantum 
Hall effect. 

In a 2D anomalous quantum Hal l'^i'^ system, the topo- 



logical edge states connect the energetically separated con- 
tinuum of energy band, and only the unidirectional topologi- 
cal edge state contribute to the conductance in the band gap. 
Due to the topological stability of Chern numbers'"* carried 
by the extended states, the conductance remains quantized in 
the presence of weak disorder. At strong disorders, the mo- 
bility edges is crossed, and MIT— occurs. However, outside 
of the gap, the bulk state and edge state may co-exist. Be- 
cause the edge state is robust for disorders, it could happen 
that the bulk state is completely killed before the edge state is 
killed, which may lead to the quantized conductance plateau. 
The quantized value is determined by the number of the robust 
edge state. Taking into account of edge state tied to both the 
spin up and spin down, disorder induced TAI can be formed, 
which is confirmed by the following calculation on the mod- 
ified Dirac model. So, roughly speaking, as long as the bulk 
states and the robust enough edge state coexist in the system, 
disorder induced quantized conductance can emerge. 

In this paper, in order to effectively study the nature of 
TAI phenomena, in addition to HgTe/CdTe quantum wellr^ 
we first find two other different models in which TAI phe- 
nomena also exists as in RefS The first model is a modi- 
fied Dirac model on a square lattice with the quadratic cor- 
rections k-cr- included and electron-hole symmetry kept. This 
model is similar to HgTe/CdTe quantum well model except 
the e-h symmetry is broken in HgTe/CdTe quantum wellji^ 
The second model is a graphene model on a honeycomb lat- 
tice with the next nearest-neighbor coupling and a staggered 
sublattice potential^. In this model the e-h symmetry and in- 
version symmetry are all broken, and the anti-directional topo- 
logical edge states tied to two opposite edges are asymmetri- 
cally distributed. Our calculation shows that in the first and 
second models, moderate disorder induces a transition from 
an ordinary metallic state to a TAI with the quantized con- 
ductance of Go = e^/Zi. The wider the ribbon is, the smaller 
the fluctuation. In the disorder induced quantized conduc- 
tance regime, the localization length of the ribbon structure 
is extremely long comparing to that of the cylinder structure. 
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FIG. 1: (Color online) A schematic diagram of the infinite long 
ribbon in square lattice for the first model [panel (a)] and the infinite 
long ribbon in honeycomb lattice along the zigzag direction for the 
second and third model [panel (b)]. The left lead and right lead (the 
red area) in the p-region [green lattice region]. 



which clearly reveals the topological origin of the transition. 
When the disorder is very strong, mobility edge crosses and 
MIT occurs. This process is related to the bulk state. In or- 
der to vividly show how the transport electron is scattered by 
the disorder, we calculate the individual local current den- 
sity from all channels including the edge state channel and 
bulk state channel. In addition to the first and second model, 
we also calculate local current density in a third model: the 
honeycomb graphene system— in which we consider Rashba 
spin-orbit interaction, exchange energy and staggered sublat- 
tice potential. In the third model, the disorder can't quantize 
the conductance but just flatten it, which is different from the 
first and second model. However, for all three models, we find 
that disorder kill all the bulk current and induce edge current 
with uni-direction. It means that in the process of transport, 
due to the topological nature, the edge state is maintained all 
the time, no matter whether it is kept in the edge channel or 
scattered into the bulk channels. So, in the system where the 
edge state and bulk state coexist, in the moderate disorder, the 
edge state is maintained while the bulk states are killed, which 
leads to the flattened conductance. If the edge state is robust 
enough to resist the disorder that is so strong that the bulk state 
are completely killed, quantized conductance can be formed. 



The rest of the paper is organized as follows. In Sec. II, 
the Hamiltonian of three model systems in the tight-binding 
representation are introduced. The formalisms for calculating 
the conductance and the local current density vector are then 
derived. Sec. Ill gives numerical results along with some 
discussions. Finally, a brief summary is presented in Sec. IV. 



II. MODELS AND FORMALISM 

A. three model Hamiltonian 

The first model is the modified Dirac model with a 
quadratic corrections k~o-,, which has the form 

Hi = ^[//^(k) + //i(k)], H^(k) = //;(-k) 

k 

//-r(k) = A(k.,o-, - kyo-y) + {m + Bk- k)o-, + £{r)o-Q (1) 

where cr^,,.,; are Pauli matrices presenting the pseudospin 
formed by s, p orbitals. In Eq.([T]), the momentum k is a good 
quantum number for periodic systems. This model is similar 
to the low-energy effective Hamiltonian of a HgTe/CrTe quan- 
tum welliS. except that the e-h symmetry is kept here while 
it is broken in HgTe/CrTe quantum well. Using substitution 
kj — > and k,, — > -idy Eq.([TJ can be transformed into the 
real space representation^^ In the calculation we consider the 
infinite long ribbon with finite width, hence the tight-binding 
representation is convenient. The tight-binding Hamiltonian 
in square lattice is given by;^ 

H-^ ^ d! (fio-o + Co--) di 

i 

+ ^ d'.^ ^fo", -H '"CTy j d\+s,y + h.c. (2) 

i 

where ctq is a unitary 2x2 matrix, e\ is a random on-site po- 
tential which is uniformly distributed in the region [w/2, w/2]. 
i = (ixJv) is the index of the discrete site of the system 
sketched in the Fig. 1(a) on the square lattice, a.v = [fl,0] and 
a, = [0, a] are the unit vectors of the square lattice with the 
lattice Constanta. d\ - [c/jj, c/pj]^ with T' denoting transpose 
and ds(p\\ and ^f j are the annihilation and creation operators 

for s(p) orbital at site i. Here C - m- 4t and t = B/a^ denote 
the nearest neighbor coupling strength. e\ in the first term of 
Eq.(|2l) is the on-site random energy account for the disorder 
The second term in Eq.© is the corrected linear Dirac term, 
in which the cr- involved terms are the quadratic corrections 
to the Dirac Hamiltonian. The individual spin up Hamiltonian 
and spin down Hamiltonian Hi in Eq.© are time rever- 
sal symmetric to each other Since they are decoupled, we can 
deal with them individually. So we shall focus only on spin-up 
Hamiltonian //| in the following calculation. 

The second model is proposed by Haldane that consid- 
ers the graphene with next-nearest neighbor coupling and 
staggered sublattice potential whose Hamiltonian can be ex- 
pressed asi^ 

Hii^) ~ [fCOc^o + toO-z + cos(k ■ ai)a-x - sin(k ■ a,)cr,,j 



+2\t„\ 



cos ^ cos(k ■ b,)cro - sin cf> ^ sin(k ■ b/)cr; 



(3) 
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where cr , . are Pauli matrix denoting the pseudospin formed 
by AB sublattice, three nearest neighbor unit vectors a,= 1,2,3 
and next-nearest neighbor unit vectors b,=i,2,3 are given by 
ai = fl(0, 1), 32 = fl(- V3/2,-l/2), V3/2,-l/2), 
bi = b(l,0), b2 = b(-l/2, V3/2), b3 = M-l/2,-V3/2) 
with a = 0.142«m and b - y/3a denoting the distance 
between nearest neighbor sites (lattice constant) and next- 
nearest neighbor sites respectively. In Eq.(l3]l, to is the stag- 
gered sublattice potential, t = 2tiVf/3a is the nearest neighbor 
coupling strength with Fermi velocity Vf - 0.89 x lO^mi"', 
the next-nearest neighbor coupling t„ - \t„\e"^ where \t„\ and 
(p are the coupling strength and phase deduced from the ef- 
fective internal magnetic field B{r) along z direction. Here f„ 
eliminates the e-h symmetry of the energy bands as shown in 
Fig. 2(b) and B{r) breaks the time-reversal invariance. Eq.([3]l 
is similar to Eq.([TJ, which includes all the cr l y,, terms that is 
A:-dependent. However, they have different symmetries, e.g., 
in Eq.lO both time-reversal and e-h symmetries are broken. 
In the tight-binding representation, Eq.([3]l can be expressed 
as; 

H2 — c/.' (fpcr. + ei)di + tc y^/flj ^i+a, + h.c.) 
i i./ 

+ if„i 2 [^'^(fl^fli+b, - /^;/?i+b,)] + (4) 

where i = (i i, i, ) is the index of the discrete site of the hon- 
eycomb lattice which includes two sublattice A (open cir- 
cle) and B (filled circle) as sketched in the Fig. 1(b). e-, - 
diag(ei, ei+ai) denotes the random potential induced by disor- 
der, d\ - [fli,fei+a,]^, fli(^i) and a!^(b'^) are the annihilation and 
creation operators for sublattice A(B) at site i. 

In the third model, we consider a graphene sheet with 
Rashba spin-orbit interaction Ae^-ikxs) as well as an exchange 
fieldJ^ In the tight-binding representation on the honeycomb 
lattice, corresponding to three directions of translational sym- 
metry, we have ki=i,2,3 —idi. Projecting them to the x andy 
direction, we can get k j and k,,. Then, the Hamiltonian of the 
third graphene model including Rashba spin-orbit coupling, 
the staggered sublattice potential and the exchange term can 
be expressed in the tight-binding representation in the follow- 
ing form>ii 

i/3 — ^ djitocr- <S> So + ffCTo ® + Ei)d; 

i 

+ ^a-[tcSo + itre^-(sxSi)]bi+^.+h.c. (5) 

ij 

where sq is the unitary matrix in spin space, s = (s^, Sy), 
Sx,y,z are Pauli matrix denoting the real spin, d^ - [ai , b'^^^^]^, 

fl| = [fljj,fl|j], bl - [^jp where fl|(|)j and ^jdjj are the 
creation operators for sublattice A and B for spin up(down) 
at site i. s; = diag(eiSo, ei+aiSo) is the random potential in- 
duced by disorder Similar to the Hamiltonian of the second 
model, tc is the nearest neighbor coupling and f,, describes the 
staggered sublattice potential. In Eq.©, t^, is the exchange 
energy which can be achieved by either a magnetic insulator 



substrate or adsorbing transition metal atoms (e.g., iron, cop- 
per) on graphene.^^ Finally tr is the strength of Rashba spin- 
orbit coupling that has been shown to be fairly strong 

B. differential current density and conductance 

In the following, we consider two geometries; ribbon and 
cylinder Fig. 1(a) and (b), respectively, depict the ribbon ge- 
ometry on square lattice and honeycomb lattice. In these sys- 
tems the finite scattering region [green (or gray in print) re- 
gion], in which random disorder is considered, is connected 
to the external reservoir through semi-infinite lead. In the fol- 
lowing, we will derive the conductance and the current density 
in the scattering region. 

For the general Hamiltonian H - Tja-i.ir'j c^Hcri,cr'jCo-'j, lo- 
cal current flowing from site i with real spin or pseudospin cr 
can be expressed asi^^ 

J<,,i(f) = -e<A?<,,i(,)) = |<c;^/f)c<.,i(f),H) 
o-'i 

- ^JiTi,tr'j (6) 
o-'j 

where e is the electron charge, G*^ ,. = i(c\.c, ,) is the matrix 
element of the lesser Green's function of the scattering region. 
Here Jo-i.o-'j is the current from site i to j. Under dc bias, the 
current is time-independent. After taking Fourier transform, 
current Jcri,o-'j can be written as; 

2e r dE r 1 
^TJ 2^^^ [Gti,^'i(E)flcr'iM] (7) 

Due to the current conservation, J^r,] = at each site in- 
side the scattering region. From now on, we will calculate the 
current density from Ja-i,a-'j- For the square lattice, it is easy 
to calculate the current density by summing over al l proj ec- 
tions of J(ji,o-'j along x and y directions as done in Ref l9l20i2l[ 
However, it is much more complicated for graphene because 
in the tight-binding representation on graphene, the current 
can flow from site i to its nearest and next nearest neighbor 
sites j. Hence we will use another simple definition?- of cur- 
rent density j - pv where p is the charge density and v is the 
velocity given by v = -i/fi{rH - Hr). The current density is 

^■^'y^^-'lj f R^Z[G'^('--'/.vH-Hr,,)]^.,^, (8) 

where r is a diagonal matrix in the discrete real space. It 
should be noted that Eq.® is valid when site i is not on the 
interface between the lead and the scattering region. For these 
boundary sites, the current density perpendicular to the inter- 
face obtained from Eq.® should be multiplied by two. This 
is due to the following reason. When we have a finite scat- 
tering, the current is conserved if the self-energy is taken into 
account. However, for these boundary sites the current from 
the lead is not accounted for by Eq.®. 
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From the Keldysh equation, the lesser Green's function is 
related to the retarded and advanced Green's functions. 



G^{E) = G'XE) 



G"iE) 



(9) 



where the sum index a - L,R denote the left and right semi- 
infinite lead, ^^(ii) = iTa{E)fa{E) in Eq.(|9|l is the lesser self 
energy of the lead-a with /„(£) = fo{E — eVa) the Fermi dis- 
tribution function. Here r„(£') = - S° ) and ll^' is related 
to the surface Green's function which can be calculated us- 
ing a transfer matrix methodi^ Va is the external bias in the 
terminal-a. In general, G^{E) can be divided into equilibrium 
and non-equilibrium parts,— 



G^(£) = G'iE) 
+ G'-{E) 



iME)Y,W) 



G"(E) 



G"(£)(10) 



where the equilibrium term can only generate persistent 
current and does not contribute to the transport<2^ so it can be 
dropped out from now on. It is the non-equilibrium term that 
gives the response to the electron injection from the source 
lead. Setting source bias = -V and drain bias Vd = 0, we 
have 



G^{E) = iG'-{E) [ME) - ME)] r,G"(£) 



(11) 



where Fj is the linewidth function of source lead. Substituting 
Eq. (fTTT l into Eq.®, the differential local current density vec- 
tor dixiyldV at site i can be expressed in the following form;^ 



(12) 



with 



p = —G^Ef + eV)Ts{EF + eV)G''{EF + eV) (13) 
2n 

Vct/v = -^(rx/vH-Hr,./,,) 
n 

where p and \c,.\/y density matrix with incident energy 
Ef + eV and velocity matrix in the central scattering region, 
respectively and we have assumed the temperature is zero. 
XXX 

When the electron is in the eigenmodes of the semi-infinite 
lead there is no scattering, the the linewidth function of 
source lead is related to the incident velocity in the form 
^Vs - U 'TsU - F^S where V j is diagonal matrix composed by 
the nonzero velocity in propagating mode and zero velocity in 
evanescent mode, U is ranked by eigenmodes including prop- 
agating mode and evanescent mode, here U can be considered 
as a unitary transformation matrix which transforms the gen- 
eral Hilbert space into the eigen channel space of lead. Then, 
f 5 can be regarded as incoming velocity matrix in source lead 
with relation v.s - (l//i)f 5 [for detail please refer to Ref lZTll . 



In the calculation, we can peak up only the propagating mode 
to construct the effective linewidth function f 5 as 



(14) 



where the sum is taken over propagating modes. A„ is the n-th 
column of matrix [U ' which is related to the «-th propagat- 
ing mode. Then from Eq.(fT2li.(fT3b, we can write the differ- 
ential local current density vector in «-th propagating eigen- 
channel: 

dl,,,/y(i)/dV ^\Yj V^^-^^''^!y\^M ^^^^ 



where 



A, = — G'-A,,„v,,„Al„G'' 



(16) 



Concerning for the scalar current flowing into the drain lead 
Jd and the conductance G - dJd/dV, we can replace \c,x/y 
with V(/ = {l/fi)td in Eq.dTSt where F^ is the effective outgo- 
ing velocity matrix of drain lead similar as effective incident 
velocity matrix v, = {l/h)Ts in source lead. Then, consid- 
ering the representation transformation, from Eq.(fT2b. ( fT3] ). 
we can get the Landauer-Buttikker formul^, which leads to 
conductance in zero temprature 



G(Ef + eV) = ^Tr [Re(t/pt/"ff rf)] ^ jT (17) 



where T = Tr [Re(G''F5G"F,/)] is the transmission coefficient 
from source lead to drain lead. 



III. NUMERICAL RESULTS AND DISCUSSION 

In the numerical calculations, the energy is measured in the 
unit of the nearest neighbor coupling constant t. For the first 
model, t - Bjcr with the square lattice constant a - 5nm. 
For the second and third model in the honeycomb lattice, 
t - 2Svir/(3a) with the carbon-carbon distance a = 0.142«m 
and the Fermi velocity \'f = 0.89 x lO^mi"' as in a real 
graphene sample.^aS, The size of the scattering region [the green 
(or gray in print) region] MxN is described by integers M and 

corresponding to the width and length, respectively. For ex- 
ample, in Fig.l, the width W - Ma with M = 4, the length 
L-Na with M = 7 in the panel (a), and the width W - Mx3a 
with M -3, the length L - N x y/3b with M = 7 in the panel 
(b). 

In all three models Eq.(l2]i, (|4]i and (|5]l, edge state exists in 
the absence of external magnetic field, which induce the quan- 
tum anomalous Hall effect. In Fig.2, we plot the band struc- 
ture (for calculation see Reflsol) and indicate the Fermi energy 
(the gray lines) in three models. We can see the unidirectional 
edge states along the each edge, down edge [the state in the 
blue (or dark gray in print) dashed lines] or up edge [the state 
in the red (or light gray in print) lines] in these models. For 
the first model [panel(a)], since e-h symmetry is kept so that 
we can focus only on electrons, i.e., only the positive energy 
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k^aW k^b(7t) k^b(7i) 



FIG. 2: (Color online) The band structure for first [panel(a)], sec- 
ond [panel(b)] and third [panel(c)] model. The horizontal gray lines 
marked the Fermi energy Ep = 0.7, 0.15 and 0.1 in three models, re- 
spectively. The blue (dark gray in print) dashed lines denote the edge 
state located in the down edge, and the red (light gray) lines is for the 
edge state located in the up edge. 



is considered. This is different from the HgTe/CdTe quan- 
tum well, in which e-h symmetry is broken and the edge state 
is more localized along the edge for the positive energy than 
the negative energy. In the calculation we set Fermi energy 
Ef = 0.7f, A/2a = 1.35f, C = 3.65f, B/a^ = t = 27.5meV, 
a = 5nm. For the second model [panel (b)], the following 
parameters are used: Fermi energy Ef - 0.15f, the nearest 
neighbor coupling constant t = 21.5eV that is used as the en- 
ergy unit in the graphene models, the next-nearest neighbor 
coupling strength t„ = O.lf * e"'^^, the staggered sublattice po- 
tential to - 0.2t. In this model, e-h symmetry is broken by 
the next-nearest neighbor coupling t„, and the edge states are 
favored the hole system (negative energy). In the third model, 
different from the second one, the next-nearest neighbor cou- 
pling is absent, the Rashba SOC = 0.1 8f and exchange en- 
ergy tg - 0.2t are considered, Fermi energy Ef = O.lf. Due 
to the Rashba SOC, the spin degeneracy is lifted and there 
are two unidirectional edge states [red or blue lines in panel 
(c)] corresponding to different spins. In the second and third 
models, we find that the staggered sublattice potential t„ is im- 
portant to study the nature of TAI phenomena. First, to breaks 
the inversion symmetry, the left [red line in Fig. 2(b), blue line 
in Fig. 2(c)] and right flowing edge states as well as the bulk 
states are asymmetrically distributed as shown in Fig. 2(b) and 
(c). As a result, the left and right injected current density are 
distributed differently, although they contribute to the same to- 
tal current. Second, to makes the bulk gap smaller, which leads 
to the coexistence of the robust edge state and bulk states in 
certain energy window AE. For instance AE - [0.1,0.5] in 
Fig. 2(b) where only one edge state near the upper edge that 
coexist with the bulk states. In Fig. 2(c), the bulk gap is closed 
completely, as will be discussed later that for this case it is 
harder to induce the TAI phenomena. However this model 
can give us some insight about TAI phenomena. 

In Fig. 3(a) and Fig. 4(a), we plot conductance and its fluc- 
tuation vs disorder for the first model in a square lattice and 
the second model in a honeycomb lattice, respectively. In 
Fig. 3(a), Fig.4(a) and Fig. 5, all data are obtained by averag- 
ing over 5000 configurations. In Fig. 3(b) and Fig. 4(b), we 
also plot the renormalized localization length ^m/W of a 2D 




disorder (t=27.5mev) 

FIG. 3: (Color online) Conductance, fluctuation and localization 
length vs disorder for the first model. Panel (a): conductance [black 
lines] and fluctuation [green or gray in print lines] in a square lattice 
with different widths W = M x a. Panel (b): renormalized local- 
ization length of a 2D strip sample [black line] and cylinder sample 
[green or gray in print line], the width of strip or the diameter of 
cylinder are all 50a. 




disorder (t=2.75eV) 

FIG. 4: (Color online) The same as Fig. 3 except that it is for the 
second model. 



strip with the width W - Ma for the square lattice [Fig. 3(b)] 
and W = Mx3a for the honeycomb lattice [Fig. 4(b)]. Further- 
more, in order to highlight the topological nature, we also plot 
the renormalized localization length ^m/W in a 2D cylinder 
sample (where the edge state is removed) with diameter W. 
Here the Fermi energy Ef is set in the bulk state [see Fig. 2]. 
In the clean system with w = 0, the conductance is contributed 
by bulk states and one edge state. When disorder is intro- 
duced, the following observations are in order: (1) the small 
disorder rapidly suppresses the conductance and enhances its 
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fluctuation. This means that the system is leading to the dif- 
fusive regime, in which the bulk extended state is scattered by 
impurities and gradually becomes localized by disorder. As a 
result, the localization length decreases accordingly for both 
strip and cylinder samples. (2) at moderate disorders, in the 
strip sample, the conductance stops decreasing and develops 
into a quantized plateau, a new phase occurs. At the same 
time, the fluctuation is deduced to almost zero, indicating the 
formation of either a bulk insulator or metallic phase. The 
abruptly increased localization length (black lines) in panel 
(b) clearly indicates the metallic phase induced by edge state. 
For the cylinder sample, however, because of the absence of 
the edge state, the system can only develop into the bulk in- 
sulator, so the localization length continues to decrease. Our 
results show that the wider the ribbon is, the larger ^m/W in 
ribbon geometry and the smaller ^m/W in cylinder geometry 
[see Fig.4(b)]. Therefore, concerning ^m/W, the peak struc- 
ture in strip geometry and the valley structure in cylinder ge- 
ometry clearly indicate the topological nature of the TAX phe- 
nomena. (3) when the disorder is strong, the quantized con- 
ductance plateau in the strip sample starts to deteriorate and 
its fluctuation begins to increase. This shows that the edge 
state will be destroyed at the strong disorder and the system 
will enter the insulating regime. 

For the cylinder sample, there is only one system trans- 
forms from band insulator to Anderson insulator. The peak 
of ^m/W indicates the transition points. (4) at the critical dis- 
order strength Wc MIT occurs. Wc is marked in Fig. (3) and (4) 
by the vertical green gray in print) dotted lines, where con- 
ductance G = 0.5^ and ^m/W ^. For cylinder sample, ^m^oo 
diverges at Wc- It should be noted that despite of the different 
topological nature, the MIT occurs at nearly the same Wc for 
strip and cylinder samples. (5) in the quantized conductance 
regime, increasing the sample size, the quantized plateau is 
much wider with smaller fluctuation. This is not surprising 
since for the larger sample, the bulk state is scattered more 
frequently while the overlap of opposite unidirectional edge 
states is smaller. This size effect is clearly seen in Fig. 5, in 
which the quantized conductance and its fluctuation in the 
midpoint [w = 9 for panel (a) and w = 2.25 for panel(b)] 
of conductance plateaus are plotted vs width of square sample 
for the fist model and second models. We can see with the 
increasing of the sample size, the conductance eventually sat- 
urates to the quantized value, its fluctuation, however, contin- 
ues to drop. It is expected that when the size is large enough, 
the complete quantized plateau will be formed with no con- 
ductance fluctuation. 

Up to now, we have confirmed that in addition to the 
HgTe/CdTe quantum well, the TAI, i.e., disorder induce quan- 
tized conductance can also occur in other systems such as the 
Dirac model and the graphene system. XXX Since the phe- 
nomenon of TAI is not the particular properties of HgTe/CdTe 
quantum well, there must be the common nature of TAI phe- 
nomena in different systems. For this reason, we study the 
scattering process in different systems through monitoring the 
local current density. 

First we study the Dirac model. Due to the inversion sym- 
metry in the Dirac model, the left and right injected current 
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FIG. 5; (Color online) Conductance [black lines] and fluctuation 
[green or gray in print lines] of a 2D square sample vs sample size in 
the first model for the fixed disorder strength w = 9 [panel (a)] and 
in the second model for the fixed disorder strength w = 2.25 [panel 
(b)]. 



are equivalent, here we focus only on the left injected cur- 
rent. We have set Fermi energy to be Ef = 0.7 for Fig. 6. In 
Fig. 6, Fig|8] Fig|8] and FigHO] all data are obtained by av- 
eraging over 1000 configurations. For Ef = 0.7, there are 
several transmission channels and we have studied the first 
three transmission channels, one is edge state and the other 
two are bulk states. In Fig.6(a)-(c), we plot the left injected 
local differential current density distribution from the lowest 
three transmission-channel including the edge channel [panel 
(c)] and the first [panel (b)] and second [panel (a)] bulk chan- 
nel in a square sample [between two vertical gray lines]. For 
the clean system with w = [the left column], the system 
is ballistic and there is no scattering. So, the current density 
from all transmission channels are distributed along the trans- 
port direction. For the left column of Fig. 6, we can see the 
right going edge state on the upper edge and the extended bulk 
state with one or two transverse peaks in the current density. 
Note that the coming eigen-channels are classified according 
to the lead. When disorder is present, due to the mode mixing 
the electron in edge channel can be scattered into bulk chan- 
nels and vice versa. As a result, the eigen-channels of the lead 
are no longer the eigen-channels of the system. Correspond- 
ing to the left column of Fig. 6, in right column, we plot the 
current density distribution at the entrance of the TAI (quan- 
tized conductance) with vv = 6. We can see that the right going 
edge current in all transmission-channels of the lead are well 
protected. 

In order to study the nature of the edge state in the pres- 
ence of disorders, we have calculated the eigen-spectrum of 
the transmission matrix FiG'^rRG".— In Fig|7jb), we plot all 
the transmissions of eigen channels for the the first and sec- 
ond model. For comparing, in Fig. 7(a), we also plot the eigen 
transmissions of HgTe/CdTe quantum well in which except 
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FIG. 6: (Color online) Left injected local differential current den- 
sity vector distribution for the lowest three sub-channels in the first 
model. The direction and the length of arrow in every site show 
the direction and the magnitude of current density vector in this site. 
Concerning a visual effect, the magnitude of current vector is also ex- 
pressed by different color. With the increasing magnitude, the color 
changes from blue to red. The left and right column are correspond- 
ing to the clear sample and dirty sample with disorder strength w = 6. 



for the terms in the first model, the term of Dk ■ k are also 
included. Here, we set D/a- - 0.75f with t = B/a^ and the 
other parameters are same as in the first model. In Fig|7] all 
data are obtained by averaging over 100 configurations. We 
found that there is always one eigen-value that is nearly equal 
to one which is identified as edge state by examining its cur- 
rent density profile. It illuminate that due to the topological 
nature, the edge state is always protected during the transport, 
although it may be scattering into the bulk channels of lead. 

From Figlll we also found that the eigen transmissions of 
the first and second model in this paper exhibit the same prop- 
erties as that in the HgTe/CdTe quantum well. It strongly im- 
plies they have the same 

Next, we examine the graphene model. We set Ef - 0.15 
for FiglS] and Fig|9] In the graphene model, the staggered 
sublattice potential is considered, which breaks the inversion 
symmetry, and the configuration of left and right injected cur- 
rent density are different. In a clean system, the left (right) 
injected current is contributed by the states that have a posi- 
tive (negative) group velocity, i.e., dE/dk^ > 0(< 0). For the 
fixed Fermi energy, the crossing of Fermi energy and energy 
band determines the momentum k of all propagating states. 
From the band structure of the graphene model [Fig. 2(b) and 
(c)], we can see that the edge state is absent for the right prop- 
agating band in the second model while in the third model, it 
is absent for the left propagating state. 

In Figl8]and Fig|9] we plot the left and right injected cur- 
rent density in the second model, respectively. We can see that 
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FIG. 7: (Color online) transmissions of system eigen channels in a 
square sample for the first model [panel (b), the width W = 100a], 
the second model [panel (c), the width W = 92* 3a], the third model 
[panel (d), the width W = 100 * 3a] and HgTe/CdTe quantum well 
[panel (a), the width W = 100a]. Different lines long the red arrow 
correspond to the eigen edge channel, the first, second, eigen bulk 
channel. 
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FIG. 8: (Color online) Panel (a-b): Left injected local differential 
current density vector distribution for the lowest two sub-channels in 
the second model. The left and right column are corresponding to 
the clear sample and dirty sample with disorder strength w = 1.8. 



for the left injected current, although edge current is absent in 
clean limit, it emerges when disorder is turned on. This can be 
understood as foIIowsjiS When the disorder is present, the ef- 
fective medium theory similar to the one in Reffiol shows that 
the parameters in the second model such as f„ and Ef are all 
renormalized and depend on disorder strength. After renor- 
malization due to the disorder, the new Fermi level (renormal- 
ized) is inside the bulk gap of the new band structure giving 
rise to a true edge state. 

Besides, the random disorder also eliminates the difference 
between states with different momentum k due to its random 
distribution. So, although the profiles of left and right injected 
current density are different in the clean system, when the 
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FIG. 9: (Color online) Panel (a-c): Right injected local differential 
current density vector distribution for the lowest three sub-channels 
in the second model. The left and right column are corresponding to 
the clear sample and dirty sample with disorder strength w = 1.8. 



FIG. 10: (Color online) Panel (a-c): Left injected local differential 
current density vector distribution for the lowest three sub-channels 
in the third model. The left and right column are corresponding to 
the clear sample and dirty sample with disorder strength w = 1.8. 



moderate disorder is considered they tend to be the same since 
the bulk state are killed and only edge state survives [see Fig.7 
and Fig. 8]. Of course, the left going and right going edge 
states are along the opposite edge. Similar to the first model, 
in clean limit, the current density in all eigenchannels are kept 
and uniformly distributed along the transport direction, except 
they have different profiles for left and right injected current 
density. 
XXX 

Form Fig. 6, [8] and |9] it seems that the phenomena of disor- 
der induce quantized conductance can appears in any system 
where quantum anomalous Hall effect is hold. However, it 
shoulder be noted that in the first and second model, the edge 
state is robust so that they can penetrate into the scattering re- 
gion under the interface scattering. In this way, the quantized 
is formed. In the following we will study another quantum 
anomalous Hall system, it is the third model. Similar to the 
second model, due to the inversion asymmetry, the configura- 
tion of the left and right injected current density are different 
in the clean limit, here we focus only on the left injected case 
in which the edge current is present in clean system. In the 
third model, the bulk energy gap is closed by the staggered 
sublattice potential, there is no pure edge state for any inci- 
dent energy, so the edge state is not so robust as in the first 
and second model. Furthermore, due to inversion asymmetry, 
one of the two unidirectional edge states turns into bulk state 
[see Fig. 2(c)]. So the disorder induced conductance must be 
e^/h if any. 

As a result, because the edge state is not so robust, the back 
following edge current is strong, although there is one unam- 
biguous unidirectional edge state [see left column of Fig fTOl . 
it is drastically reflected by the interface scattering, then the 



edge current can't deeply penetrate into the scattering region 
when moderate disorder strength is considered, which can be 
clearly seen in FigfT0l'c2). On the other hand, the edge state 
that enter into the scattering region can always survive in a 
moderate disorder [see Fig fT0l" a2) and (b2)]. As a result, the 
conductance can't be quantized but flattened by the moderate 
disorder, which is different from the first and second model, 
but they share the same nature, that can also be seen from 
Figlltc) and (d). 

In FiglTJc) and (d) we plot the system eigen transmissions 
of the second and third model. We can see for all model (the 
first, second, third model and the HgTe/CdTe quantum well), 
the transmissions of eigen bulk state are rapidly reduced to 
zero in a weak disorder, while the eigen edge state can be 
kept in the moderate disorder. When the eigen edge state is 
robust, it can be completely kept [transmission equal to one, 
see panel (a), (b) and (c)], but when the eigen edge state is not 
so robust, it can only be partly kept [transmission is smaller 
than one, see panel (d)]. So we can conclude that in the system 
where robust edge state and bulk state coexist, the edge state 
is maintained while the bulk states are killed in the moderate 
disorder, which leads to the flat conductance. If the edge state 
is robust enough to resist the disorder that is so strong that the 
bulk state are completely killed, quantized conductance can 
be formed and TAI phenomena appears. 



IV. CONCLUSION 

In summary, the disorder effect is studied in several quan- 
tum anomalous Hall system. When beyond the energy gap. 
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the edge state and bulk state often coexist in clean system 
or weak disorder system. In a moderate disorder, the inter- 
band scattering happens. If the edge state is robust enough 
to resist the interface scattering, it can penetrate deep into 
scattering region where the edge state and bulk state coexist 
in all of the coming eigen channels due to the mod mixing. 
Due to the topological nature, the bulk states are drastically 
killed and edge state survives. In the eigen spectrum of the 
whole system, the eigen channel corresponding to the edge 
state is kept well and the other bulk channels are completely 
killed, which leads to a flattened conductance accompanied 
with the rapidly decreased fluctuation. In several spatial sys- 
tems such as HgTe/CdTe quantum well, modified Dirac model 
and the graphene system with next-nearest coupling and stag- 



gered sublattice potential, the edge state is very robust so that 
it can survive in a disorder that is so strong that the bulk states 
are completely killed by them. Without the bulk states, the 
unidirectional edge state located in up or down edge can't be 
transformed to the opposite edge, the the back following cur- 
rent is then prohibited. As a result, the quantized conductance 
can be formed and TAI phenomena that disorder induce quan- 
tized conductance appears. 
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